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ABSTRACT 

We use the public release of ~ 100, 000 galaxies of the Two Degree Field Galaxy 
Redshift Survey (2dFGRS) to analyze the internal dynamics of galaxy clusters. We 
select 43 non-interacting clusters which are adequately sampled in the 2dFGRS public 
release. Members of these clusters are selected out to ~ 2 virial radii. We build 
an ensemble cluster by stacking together the 43 clusters, after appropriate scaling of 
their galaxy velocities and clustercentric distances. We solve the Jeans equation for the 
hydrostatic equilibrium for the member galaxies within the virial radius of the ensemble 
cluster, assuming isotropic orbits. We constrain the cluster mass profile within the virial 
radius by exploring parameterized models for the cluster mass-density profile. We find 
that both cuspy profiles and profiles with a core are acceptable. In particular, the 
concentration parameter of the best fit NFW model is as predicted from numerical 
simulations in a ACDM cosmology. Density profiles with very large core-radii are ruled 
out. Beyond the virial radius, dynamical equilibrium cannot be taken for granted, and 
the Jeans equation may not be applicable. In order to extend our dynamical analysis 
out to ~ 2 virial radii, we rely upon the method which uses the amplitude of caustics in 
the space of galaxy clustercentric distances and velocities. We find very good agreement 
between the mass profile determined with the caustic method and the extrapolation to 
~ 2 virial radii of the best-fit mass profile determined by the Jeans analysis in the 
virialized inner region. We determine the mass-to-number density profile, and find it is 
fully consistent with a constant within the virial radius. The mass-to-number density 
profile is however inconsistent with a constant when the full radial range from to ~ 2 
virial radii is considered, unless the sample used to determine the number density profile 
is restricted to the early- type galaxies. 

Subject headings: galaxies: clusters: general - galaxies: kinematics and dynamics - dark matter 
- cosmology: observations 



1. INTRODUCTION 

The determination of cluster masses has a 
long story (see, e.g. Biviano 2002, for a re- 
view). It dates back to Zwicky (1933, 1937)'s 
and Smith (1936) 's preliminary estimates of 
the masses of the Coma and the Virgo cluster. 
These early estimates were based on the virial 
theorem, using galaxies as unbiased tracers of 
the cluster potential (i.e. the light traces mass 
hypothesis). The & White (1986) and Mer- 
ritt (1987) were among the first to point out 
that cluster masses were in fact poorly known, 
since relaxing the light traces mass assump- 
tion widens the range of allowed mass models 
considerably. However, the limited amount 
of redshift data made it very difficult, if not 
impossible, to constrain the relative distribu- 
tions of cluster mass and cluster light. 

With the advent of multi-object spec- 
troscopy, a large number of redshifts for clus- 
ter galaxies became available. In particular, 
the two main catalogs of redshifts of clus- 
ter galaxies became available, viz. the ESO 
Nearby Abell Cluster Survey (ENACS, Kat- 
gert et al. 1996, 1998), and the Canadian Net- 
work for Observational Cosmology (CNOC, 
Yee et al. 1996; Ellingson et al. 1998). More- 
over, field galaxy surveys also contributed to 
increase the number of redshifts for cluster 
galaxies, in particular, the Century Survey 
(Wegner et al. 2001), and the Two Degree 
Field Galaxy Redshift Survey (2dFGRS, Col- 
less et al. 2001; De Propris et al. 2002, here- 
after DP02). 

These new data-bases prompted new more 
detailed investigations in the issue of the clus- 
ter mass determination. With a significant 
number of galaxy redshifts per cluster, it be- 
came possible to determine the cluster mass 
profile by solving the Jeans equation (e.g., 
Binney & Tremaine 1987) for the equilibrium 
of member galaxies in the cluster gravitational 
potential. Often, several cluster samples have 



been combined to improve the statistics, by 
making the implicit assumption of cluster ho- 
mology. 

From the analysis of ~ 1000 galaxies in 
the CNOC clusters, combined to form a sin- 
gle cluster, Carlberg et al. (1997a) concluded 
that galaxies trace the mass to within ±30%, 
and that the cluster mass profile is well de- 
scribed by a Navarro, Frenk & White (1997, 
NFW hereafter) profile, or by a Hernquist 
(1990) profile. Their results have substan- 
tially been confirmed by the more detailed 
analysis of van der Marel et al. (2000, here- 
after vdMOO). vdMOO found that the mass- 
to-light ratio is almost independent of radius, 
and the NFW profile is an adequate fit to the 
data, but other density profiles are equally ac- 
ceptable, depending on the orbital anisotropy 
of cluster galaxies. Carlberg et al. (2001) have 
recently applied the technique of Carlberg et 
al. (1997a) to a sample of ~ 700 galaxies in 
~ 200 groups from the CNOC2. They found 
a strongly increasing mass-to-light ratio with 
radius, a result at variance with that found 
by Mahdavi et al. (1999). A detailed analysis 
of the mass and anisotropy profiles of clusters 
from ENACS is ongoing (Katgert, Biviano &: 
Mazure, in preparation). 

A new technique for the determination of 
cluster mass profiles has recently been intro- 
duced by Diaferio & Geller (1997) and Diafe- 
rio (1999, hereafter D99). It is based on the 
determination of the amplitude of the caustics 
in the space of velocities and clustercentric 
radii. Applications of this technique include 
Geller, Diaferio & Kurtz (1999)'s and Rines et 
al. (2001) 's determination of the Coma clus- 
ter mass profile out to 14 h' 1 Mpc, Rines et 
al. (2000) 's determination of the mass profile 
of Abell 576 out to 6 h~ x Mpc, and Reiseneg- 
ger et al. (2000) 's determination of the mass 
of the Shapley supercluster. The NFW pro- 
file was found to provide a consistent fit to 
the mass profiles. Interestingly, Rines et al. 
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(2000) found a decreasing mass-to-light profile 
in Abell 576 while Rines et al. (2001) found a 
constant mass-to-light profile in Coma. 

From a theoretical point of view, both de- 
creasing and increasing mass-to-light profiles 
have been predicted. A decreasing mass-to- 
light profile could result from the tidal strip- 
ping of galaxy halos in the cluster centers 
(Mamon 2000). An increasing mass-to- light 
profile could instead result from the combined 
effects of dynamical friction and galaxy merg- 
ing (Fusco-Femiano & Menci 1998). In gen- 
eral, galaxies are unlikely to be distributed 
exactly like the mass, if anything because dif- 
ferent galaxy populations have different dis- 
tributions (e.g., Biviano et al. 2002). 

In this paper we determine the mass pro- 
file of an ensemble cluster built from the com- 
bined data of 43 clusters extracted from the 
2dFGRS (Colless et al. 2001), using the tech- 
niques described by vdMOO and Diaferio & 
Geller (1997, see also D99). Since the 2dF- 
GRS is a field survey, only a small fraction 
of the 2dFGRS galaxies are cluster members. 
Our final sample contains 1345 cluster mem- 
bers only, a number comparable to that of 
cluster members in the CNOC and smaller 
than that of cluster members in the ENACS. 
The main advantage of using the 2dFGRS for 
the determination of the cluster mass profile 
is the possibility of sampling the cluster dy- 
namics to a large distance from the cluster 
center. 

This paper is organized as follows. In § 2 
we describe the data sample, how we assign 
the cluster membership and how we combine 
the 43 cluster data-sets into a single ensemble 
cluster. In § 3 we determine the ensemble clus- 
ter mass profile. In § 4 we use this mass profile 
to determine the mass-to-number density pro- 
file. In § 5 we discuss our results and compare 
them with previous determinations of cluster 
mass and mass-to-light profiles. Finally, in § 6 
we present a summary of our results. A Hub- 



ble constant of 100 h km s Mpc is used 
throughout. 

2. THE DATA SAMPLE 

The data sample we use is the 2dFGRS 
public release version of June, 30^, 2001. 
It contains ~ 100, 000 redshifts and spectral 
types as well as 6j photometry for ~ 460, 000 
objects. The bj magnitudes are extinction- 
corrected total magnitudes derived from up- 
dated versions of the original APM (Auto- 
mated Plate Measuring machine) scans (see 
Colless et al. 2001; Norberg et al. 2001). The 
target galaxies for the 2dFGRS were selected 
to have bj < 19.45. Within the 2dFGRS 
public release, we only consider those galaxies 
in the fields of galaxy clusters. Galaxy clus- 
ters have been identified in the 2dFGRS by 
DP02, via a cross-correlation of the 2dFGRS 
with the cluster catalogs of Abell, Corwin, & 
Olowin (1989, hereafter, ACO), of Lumsden 
et al. (1992, the Edinburgh-Durham Cluster 
Catalogue, EDCC) and of Dalton et al. (1997, 
the APM cluster catalog). From DP02's sam- 
ple we exclude those clusters whose central 
regions are not evenly sampled in the 2dF- 
GRS public release, and those which, accord- 
ing to the authors themselves, have less than 
20 members with redshifts in the 2dFGRS. 
Of course, we also get rid of double entries 
in DP02's list (such as, e.g., APM 309 which 
corresponds to ACO 3062). We are thus left 
with 91 clusters. 

Following Carlberg et al. (1997a) we de- 
fine r u = y/3a p /[10H(z)] ~ r2oo 5 the radius 
at which the mean interior overdensity is 200 
times the critical density (in a 0,q = 1 Uni- 
verse), usually called the 'virial radius'. As a 
first estimate of the clusters projected veloc- 
ity dispersions along the line-of-sight, <t p 's, we 
use DP02's values. 

To these 91 clusters, we apply the selection 
procedure of cluster members of Fadda et al. 
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(1996) and Girardi et al. (1998). First, we ap- 
ply Pisani (1993) 's one-dimensional algorithm 
based on the adaptive kernel technique (see 
also Appendix A of Girardi et al. 1996). The 
adaptive kernel technique is a nonparamet- 
ric method for the evaluation of the density 
probability function underlying an observa- 
tional discrete data set. The method returns 
a list of all density peaks detected, as well as 
their probabilities, and the objects associated 
to each peak. We apply this method to the 
velocity distribution of all galaxies projected 
within 2r u of each cluster center. We select as 
significant those peaks with > 98% probabil- 
ity. When adjacent peaks have 20% or more 
of their objects in common, we join them to- 
gether if they are closer in velocity than 1500 
km s _1 . When one peak is closer than 1500 
km s _1 to two other peaks, themselves sepa- 
rated by more than 1500 km s _1 , we join the 
central peak to the one which is closest, and 
leave the third peak separate. 

Quite often, several significant density 
peaks are found in a single cluster region. 
Among the significant peaks, we choose that 
one which is closest to the mean cluster veloc- 
ity, as given by DP02. In this way we avoid 
possible mis-identifications, since sometimes 
the search region around each cluster (a cir- 
cle of 2r u radius) is wide enough to contain 
another cluster. 

We identify 75 clusters with a significant 
density peak in the velocity space and at 
least 20 galaxies associated to that peak. In 
one case (ACO-S1155), no significant peak is 
found, and in another 15 cases, there are less 
than 20 galaxies in the peak corresponding to 
the cluster. These 16 clusters are excluded 
from our sample. 

We re-define the centers of the 75 clusters. 
We define the new cluster center as the po- 
sition of the highest density peak found with 
a two-dimensional adaptive kernel technique. 
In order to avoid mis-identifications, we only 



use galaxies less distant than 0.5 h Mpc 
from the old cluster center listed by DP02. 
Using these new center determinations, we ap- 
ply the 'shifting-gap' method (Fadda et al. 
1996) to get rid of the remaining interlopers. 
The shifting-gap method makes use of both 
galaxy velocities and clustercentric distances. 
Galaxies are sorted in order of increasing dis- 
tance from their cluster centers. Among sam- 
ples of galaxies in (overlapping and shifting) 
bins of 0.4 h^ 1 Mpc from the cluster cen- 
ter (or wide enough to contain > 15 galaxies 
each), gaps > 1000 km s _1 are identified in 
the galaxy velocity distribution. These gaps 
define the edges of the galaxy velocity dis- 
tribution. Outliers from this distribution are 
flagged as interlopers. 

Having excluded the interlopers, we make 
a new estimate of the cluster mean velocities, 
TT's, and velocity dispersions, <7 p 's, (and hence 
r M 's). The velocity moments are estimated 
with the robust methods of Beers, Flynn, & 
Gebhardt (1990, see also Girardi et al. 1998). 
Velocity dispersions are corrected for the ve- 
locity errors and moved to the cluster rest- 
frames (Danese et al. 1980). These new v, a p 
and r u estimates are the starting values for 
an iterative procedure. We consider only the 
cluster members within a distance < r u from 
the cluster center to recalculate v and a p . If 
the difference between the new estimate of a p 
and the previous one is < 10%, the procedure 
is halted. Otherwise, we iterate, by recalcu- 
lating v and (j p on the members within a dis- 
tance of < r u from the cluster center, where 
the new value of r u is used. With this itera- 
tive procedure, we eliminate 7 of the initial 75 
clusters, because they are left with less than 
5 members within r u at any of the iteration 
steps. 

With the new center and r u determina- 
tions, we then iterate the adaptive-kernel pro- 
cedure and the shifting-gap method for the 
selection of cluster members. This time, we 
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extend the search region to 3r u , in order to 
ensure that we will be covering a region of 
at least 2r u radius when the final definitive 
estimate of r u is done. On the selected galax- 
ies we re-determine v and a p via the iterative 
procedure described above. We are left with a 
sample of 67 clusters, after eliminating one be- 
cause it contains less than 5 members within 
ru- 
in order to determine the mass profile of 
galaxy clusters to sufficiently large distances, 
it is wise to exclude interacting clusters, whose 
dynamics is difficult to model. Say vl, <Ji, and 
r Uj i are, respectively, the mean velocity, ve- 
locity dispersion, and virial radius of cluster 
i, and Rij is the projected distance between 
the centers of clusters i and j; we say cluster i 
and cluster j are 'interacting' if the following 
two conditions apply: 

| Vi-Vj |< 3((7j + (7j) (1) 

Ri,j <2(r u ,i + r u ,j). (2) 

We eliminate 24 interacting clusters, and we 
are thus left with a sample of 43 reasonably 
isolated clusters. 

Before we can construct azimuthally aver- 
aged profiles, we must verify that the circu- 
lar region of 2r u radius around each cluster 
center is fully covered in the 2dFGRS pub- 
lic release data sample. This is indeed the 
case for all the clusters in our sample, except 
two, APM 715 and ACO 892. For these two 
clusters we define an angular sector which ex- 
cludes the region not covered in the 2dFGRS 
public release data sample, and we consider 
only the galaxies within this sector. 

Using a two-dimensional Kolmogorov-Smirnov 
test (e.g. Fasano & Franceschini 1987) we also 
verify that the spatial distributions of galax- 
ies with and without redshifts are not signifi- 
cantly different, at least out to 2r M . There are 
only four clusters for which there is a marginal 
(> 90%) evidence of a difference: ACO 957, 



ACO 1750, ACO 2734, EDCC 153. For these 
clusters the redshift completeness is not ho- 
mogeneous within 2r n , and decreases outside. 
Only 9% (respectively 2%) of all the galax- 
ies in our sample occupy regions where the 
completeness in redshift is significantly below 
(respectively 50%) the average completeness 
for their cluster. We therefore deem it unnec- 
essary to correct for this incompleteness. In 
summary, thanks to the quality of the 2dF- 
GRS, our final sample is almost unaffected by 
problems of radial incompleteness. 

Our final sample contains 3602 field galax- 
ies and 1345 cluster members with redshifts, 
within 2r u from the centers of 43 clusters. The 
absolute 6j magnitudes of the 1345 cluster 
members, k-corrected as described by Madg- 
wick et al. (2002), range from M bj = -22.3 + 
51og/i to M bj = -14.4 + 51og/i. The dif- 
ferent clusters are sampled down to differ- 
ent absolute magnitude limits, since they lie 
at different redshifts, from a limiting M bj = 
-14.4 + 51og/i in ACO S301, to a limiting 
M bj = -19.3 + 5 log h in APM 294. 

The 43 clusters have an average a p of 490 
km s _1 . This is significantly lower than the 
average a p estimated by DP02 for the same 
clusters, 776 km s _1 . Interestingly, our veloc- 
ity dispersion estimates and those of DP02 are 
significantly correlated (97% probability ac- 
cording to a Spearman rank correlation test), 
which suggests that the difference in a p is sys- 
tematic. 

We believe that DP02's a p s are systemati- 
cally overestimated, for the following reasons. 
First, DP02's highest value of a p among the 
43 clusters is 1620 km s _1 , higher than any 
of the a p values found in the complete cluster 
sample of Mazure et al. (1996). For compari- 
son, the hottest X-ray cluster known, 1E0657- 
56, has a velocity dispersion of 120llg2° 
km s _1 (Barrena et al. 2002). Second, 
DP02's Cp-estimates are too large when com- 
pared to the cluster richnesses. This can be 
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seen by considering the 27 ACO clusters in 
our sample, whose average richness count is 
41. The 27 ACO clusters have an average 
o v of 493 km s _1 according to our analysis, 
or 822 km s _1 according to DP02. The for- 
mer value is as expected on the basis of the 
cluster richness- velocity dispersion correlation 
(Fadda et al. 1996), while the latter value is 
certainly too high. Finally, our a p -estimates 
are in better agreement than those of DP02, 
with the values given in the literature for five 
of the 43 clusters (Mazure et al. 1996; Girardi 
et al. 1998; Alonso et al. 1999). The mean 
difference between our <r p -estimates and those 
of the literature is ± 172 km s _1 , that be- 
tween DP02's dp-estimates and those of the 
literature is +378 ± 101 km s -1 . 

Our cluster sample is presented in Table 1. 
We list in column (1) the cluster name (as 
given in the list of DP02), in columns (2) and 
(3) the coordinates of the cluster center, in 
column (4) the number of member galaxies 
within 2r u , in column (5) the number of mem- 
ber galaxies within r u , in column (6) the mean 
cluster velocity, and in column (7) the cluster 
velocity dispersion (computed using the ro- 
bust estimators of Beers et al. 1990). 

3. THE MASS PROFILE 

In order to put significant constraints on 
the cluster mass profile we need a sufficiently 
large data-set. We therefore join together the 
43 clusters into a single ensemble cluster. To 
this aim, we scale the projected distances of 
galaxies from their cluster centers, R's, by r u 
(see § 2) which is an approximate estimate of 
7*200 • We scale the line-of-sight velocities of 
the galaxies relative to the cluster mean ve- 
locity, (v — v), with the global velocity dis- 
persions dp of their parent clusters. Similar 
scalings have often been used in the litera- 
ture (e.g. Biviano et al. 1992; Carlberg et al. 
1997a). 
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Fig. 1. — Distribution of galaxies of the en- 
semble cluster in the space of normalized clus- 
tercentric distances and velocities. Top panel: 
large dots indicate cluster members. Bot- 
tom panel: large dots indicate the 642 cluster 
members, late-type spirals excluded, within 
R < r u, on which the Jeans analysis is based. 
The vertical dashed line indicates the R = r u 
limit beyond which the cluster membership 
assignment is less reliable. 
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Fig. 2. — Projected number density profile 
of the cluster members (late-type spirals ex- 
cluded), determined through the MAPELN 
method. Error bars were determined from 
1000 bootstrap resamplings (dashed lines). 
The dash-dotted line represents the best-fit 
to the observed profile with a function of the 
form suggested by vdMOO. 
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In Figure 1 (top panel) we show the dis- 
tribution of galaxies of the ensemble clus- 
ter in the space of normalized clustercentric 
distances and velocities. Large dots indi- 
cate cluster members. It is easy to spot a 
few apparent outliers that have been classi- 
fied as cluster members, mostly at large dis- 
tances from the cluster center. These possible 
misidentifications are due to the difficulty of 
the shifting- gap procedure in rejecting inter- 
lopers in regions of low galaxy density. In or- 
der to avoid possible problems of wrong mem- 
bership assignments, we do not use the mem- 
bership assignment beyond r u in our analysis. 

3.1. The Jeans approach 

The determination of the mass profile can 
in principle be obtained by a straightforward 
application of the Jeans equation. A rea- 
sonable assumption must however be made 
on the orbital anisotropy of the tracers of 
the gravitational potential. Several analyses 
have reached the conclusion that the orbits of 
galaxies in clusters are quasi-isotropic (Carl- 
berg et al. 1997b; van der Marel et al. 2000), 
except for late-type spirals, which are found to 
be on moderately radial orbits (see, e.g., Mohr 
1996; Biviano et al. 1997, see also Adami, Bi- 
viano, & Mazure 1998). By excluding late- 
type spirals from our sample we are therefore 
confident that we can assume isotropic orbits 
for the remaining cluster galaxies. In order to 
exclude late-type spirals we impose an upper 
limit on the spectral type parameter provided 
by the 2dFGRS, r? < 1 (see Madgwick et al. 
2002). 

As explained in § 3, our membership as- 
signment may not be very robust beyond r u . 
We therefore restrict the Jeans analysis to the 
sample of cluster members with R < r u . Ad- 
ditional motivations for this choice are: 

• reducing the effect of possible anisotropies. 
In fact, galaxies infalling into a cluster 



are likely to lose their radial anisotropy 
when crossing the denser central regions 
(Mamon 1995). 

• Reducing the influence of possible sub- 
clustering. In fact, the central cluster 
regions are a hostile environment for 
the survival of substructures (Gonzales- 
Casado et al. 1994). 

• Satisfying the condition of dynamical 
equilibrium, needed for the application 
of the Jeans equation. In fact, r u w r2oo 
delimits the region of virialization. 

In total, we have 642 cluster members with 
R < r u and rj < 1. Their distribution in 
the space of clustercentric distances and ve- 
locities is shown in Figure 1 (bottom panel). 
Their normalized velocity-distribution is not 
significantly different from a Gaussian, ac- 
cording to a Kolmogorov-Smirnov test. This 
suggests that our assumption of isotropy is 
acceptable (see Merritt 1987; vdMOO). We 
solve the Jeans equation under the isotropic 
assumption for the sample of 642 cluster mem- 
bers. For this, we follow the method of van 
der Marel (1994) and vdMOO. 

First, we determine the projected number 
density profile, S(i?) of cluster members out 
to 2r u (late-type spirals excluded) with the 
Maximum Penalized Likelihood (MAPELN) 
method of Merritt & Tremblay (1994). We fit 
the observed profile with a multi-parameter 
function of the form suggested by vdMOO (see 
eq. [2] in vdMOO, and Figure 2). The function 
is sufficiently general as to allow a very close 
representation of the real S(i?). We Abel- 
invert the S(i?) best-fitting function to pro- 
duce the 3-dimensional number density pro- 
file v{r). Abel-inversion requires the knowl- 
edge of the asymptotic behaviour of S(i?) to 
large radii. We try different extrapolations of 
the observed T,(R) to 10r u and choose the one 
that gives the v{r) which most closely match 
the obsered S(i?) after Abel-projection. 
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We then consider two functional forms for 
the mass density profile. One is taken from 
vdMOO: 

p{r)=p Q {r/a)^{l + r/af-\ (3) 

The other, often employed in the literature 
(e.g. Girardi et al. 1998), is the /3-model (Cav- 
aliere & Fusco-Femiano 1978) density profile: 

p(r)=p [l + (r/r c fr^ 2 . (4) 

Hereafter we refer to these two mass-density 
models as the £- and the /3-model. 

We consider a grid of values for the mass 
density parameters (a, £) and (r c , (3). For each 
set of values we compute the shape of the pre- 
dicted cr p and the density normalization po 
that minimizes the difference between the pre- 
dicted a p and the observed (binned) a p , out 
to R = r u . We then constrain the allowed val- 
ues of the density profile parameters in a x 2 
sense. The observed cr p -profile is determined 
with the biweight estimator (Beers et al. 1990) 
in 11 radial bins, each containing 57 galaxies. 

The allowed 68% confidence level (c.l. here- 
after) isocontours in the space of (a,£), and 
(r c , (3) parameters are shown in Figure 3. The 
confidence levels are determined as described 
in Avni (1976), for two interesting param- 
eters. Clearly, density models with a wide 
range of parameters provide an acceptable fit 
to our data. In particular, the best fit to 
the £-model is obtained for a = 0.33 r u and 
£ = 1.4, with a reduced x 2 = 0-6- A stan- 
dard NFW model (corresponding to the case 
£ = 1.0) provides an almost equally good fit 
to the data. The best fit NFW model is ob- 
tained for a = 0.18 r u , as expected for halos 
with a v ~ 500 km s~ x in a ACDM Universe. 
Also acceptable are the density models with 
a core (£ = 0), in which case the scale is con- 
strained to be small, a ~ 0.1 r u . Consistently, 
the allowed values of r c for the /3-model are 
also rather small. The best-fit is obtained for 
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Fig. 3. — Isocontours of allowed parameter 
values for the mass density £ and (5 models - 
see text. Crosses indicate the best-fit values. 
Top panel: x 2 68% c.l. isocontours in the 
space of (a, £) parameters. Bottom panel: x 2 
68% c.l. isocontours in the space of (r c , /3) 
parameters. 
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Fig. 4. — Predicted velocity dispersion pro- 
files from the best-fit mass density models, 
and the observed a v (binned points with 68% 
error bars). Solid line: £-model (a = 0.33 r u 
and £ = 1.4); dash-dotted line: NFW-model 
(a = 0.18 r u and £ = 1); dashed line: /3-model 
(r c = 0.06 r u and (3 = 0.8). 
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Fig. 5. — Mass profiles determined by the 
Jeans analysis. Solid line: mass profile ob- 
tained from the best-fit £ density model; dash- 
dotted line: mass profile obtained from the 
best-fit NFW density model; dashed line: 
mass profile obtained from the best-fit (3- 
model. Dotted lines: 68% c.l. for the mass 
profiles, derived from integration of the den- 
sity profiles which provide an acceptable fit 
(to within the 68% c.l.) to the observed a p 
(see Figure 3). 



r c = 0.06 r u and (3 = 0.8, which are close to 
the values found by Girardi et al. (1998) for 
the number density profiles of cluster galax- 
ies. Core-radii r c > 0.2 r u (r c > 0.3 r u ) are 
excluded at the 68% (90%) c.l.. A standard 
King (1962) profile ((3 = 1.0) with r c = 0.14 r u 
is also an acceptable fit; these values are close 
to those obtained by Adami et al. (1998) for 
the number density profiles of cluster galaxies. 

In Figure 4 we show the predicted veloc- 
ity dispersion profiles for the best-fit £- (solid 
line), £ = 1 NFW (dash-dotted line), and 0- 
(dashed-line) models, as well as the observed 
(binned) <7 p -profile. 

We obtain the 68% c.l. for the cluster mass 
profile by integrating the mass density pro- 
files which provide 68% c.l. acceptable fits 
to the observed <7 p -profile. The mass pro- 
files determined via integration of the best- 
fit mass density profiles, and the 68% c.l. for 
the cluster mass profile are shown in Figure 5. 
The mass is in the normalized units of our 
ensemble cluster, i.e. r u <jp/G, which corre- 
sponds to ss 5 x 10 13 h^ 1 M & , for the aver- 
age value of a v ~ 500 km ■sT 1 of the 43 clus- 
ters that constitute our ensemble cluster. We 
call M n the normalized mass. The mass pro- 
files obtained from the best-fit £- and f3- mass 
density models are very similar in the fitted 
region (r < r u ). The total mass at r < r u 
appears to be determined with an accuracy of 
~ ±8%. 

3.2. The caustic approach 

In order to extend the determination of 
the cluster mass profile beyond the virializa- 
tion radius r u ~ r2oo, we must rely on meth- 
ods that do not assume the cluster dynamical 
equilibrium, since this is unlikely to hold in 
the external cluster regions. The method we 
use is that of Diaferio & Geller (1997) and 
D99. These authors have shown that the ve- 
locity field in the cluster outskirts is deter- 
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mined by the cluster mass distribution, via 



GM{< r) - GM{< r ) = I A 2 {x) Tp(x) dx, 

Jr 

(5) 

where A is the amplitude of the caustics de- 
scribed by the galaxy distribution in the space 
of normalized clustercentric distances and ve- 
locities, and J-p is a function of the gravita- 
tional potential and the anisotropy profile. 

We determine the caustics with an adap- 
tive kernel method, as described in D99 (see 
also Pisani 1996). Following Diaferio & Geller 
(1997) and Geller et al. (1999), we scale the 
smoothing-window sizes along the axes of 
clustercentric distance and velocities before 
applying the adaptive kernel method, in such 
a way as to give equal weights to the typi- 
cal uncertainties of normalized clustercentric 
distances and velocities. In our sample the 
median velocity uncertainty is 90 km s -1 , or 
0.2 in normalized units. The positional uncer- 
tainty of a galaxy is of the order of the galaxy 
size, i.e. about 20 h^ 1 kpc, or 0.02 in normal- 
ized units. We therefore take q = 10 as the 
ratio between the smoothing-window sizes on 
the two axes. Geller et al. (1999) have shown 
that the mass profile determination is largely 
independent on the precise choice of q. As in 
Geller et al. (1999), we symmetrize the caus- 
tics with respect to the mean cluster velocity. 

Following D99 and Geller et al. (1999), we 
choose the caustic that minimizes the func- 
tion: 

S=\^ C -^{<R)\ 2 . (6) 

The v1 sc term is determined from the caus- 
tic amplitude, and f„(< R) is the normal- 
ized velocity dispersion of the cluster mem- 
bers contained within R, the average cluster- 
centric distance of all cluster members. As 
described in § 3.1 our fiducial sample of clus- 
ter members does not include late-type spirals 
and is restricted to galaxies within r u . We 
show in Figure 6 the caustic that minimizes 
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0.0 0.5 1.0 1.5 2.0 

R/r u 

Fig. 6. — 'Optimal caustic' and the distri- 
bution of galaxies from the ensemble cluster 
in the space of normalized clustercentric dis- 
tances and velocities. Large dots indicate the 
cluster members used to determine v^{< R) - 
see text. 
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0.0 0.5 1.0 1.5 2.0 

Fig. 7. — Function J-a for the determination 
of the cluster mass profile with the caustic 
method. The dashed line shows the value 
= 0.5 used by Geller et al. (1999). 
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the S function (see eq. 6). Hereafter, we refer 
to this caustic as the 'optimal caustic'. 

In order to translate the optimal caustic 
amplitude into a mass profile, we must know 
the function This function actually de- 
pends on the gravitational potential one is 
willing to determine. However, numerical sim- 
ulations show that this dependence is small 
at sufficiently large distances from the cluster 
center (D99). Diaferio & Geller (1997) and 
Geller et al. (1999) used the simple approxi- 
mation T$ = 0.5, but this is not really sup- 
ported by the results of D99 (see Figure 3 in 
D99's paper). In our analysis we adopt a non- 
constant Tr (see Figure 7) which is a smooth 
approximation of the result obtained by D99 
for a ACDM cosmology. The resulting mass 
function is shown in Figure 8 (top panel) . The 
c.l.'s on this mass function are determined as 
in Geller et al. (1999). For the sake of com- 
parison, we also show the mass function one 
would obtain by adopting T$ = 0.5. 

As can be seen from the top panel of Fig- 
ure 8, the mass profiles determined from the 
same optimal caustic, using different T$ func- 
tions, can differ considerably. As a matter of 
fact, the choice of T$ constitutes the dom- 
inant source of uncertainty in this method. 
Note however that the mass profile deter- 
mined using a non-constant is in remark- 
able consistency with the mass profile deter- 
mined via the Jeans analysis (see § 3.1). 

In order to reduce the systematic errors in 
the mass profile, it is better to avoid using 
the caustic method in the central cluster re- 
gions. In fact, the main advantage of the caus- 
tic method lies in its applicability beyond the 
virialized central region, where classical meth- 
ods fail, but classical methods may prove su- 
perior in the central region, where the caus- 
tic method suffers from the uncertainty in the 
function Tp. We can therefore take advan- 
tage of the constraints imposed by the Jeans 
analysis on the mass interior to r u , and com- 



pute the mass from the caustic method only 
at r > r u (i.e. we set r$ = r u in eq. 5). By do- 
ing this we reduce the uncertainty in the mass 
profile determined via the caustic method, be- 
cause the radial dependence of is much 
smaller at large radii than in the central re- 
gion (see Figure 2 in D99). The combined 
mass profile from the Jeans and the caustic 
approach is shown in the bottom panel of Fig- 
ure 8. Clearly, the mass profiles determined 
with different choices of Tp are now in good 
agreement, since we compute them only for 
r > r u . The uncertainty on the total mass at 
r = 2r u is « ±15%. 

We note that the mass profile determined 
by the caustic method follows very closely the 
extrapolation of the best-fit ^-profile out to 
r = 2r u . The best-fit /3-profile is somewhat 
steeper, but well within the 68% confidence 
region of the caustic mass profile. 

4. THE MASS-TO-NUMBER DEN- 
SITY PROFILE 

In § 3.2 we have determined the average 
cluster mass profile out to r = 2r u . We obtain 
the average cluster mass density profile, p, by 
numerical differentiation of the mass profile 
shown in the bottom panel of Figure 8 (in- 
dicated by the solid line up to r = r u , and 
by the dash-dotted line at larger radii). We 
determine the 68% confidence interval on this 
density profile by constructing a set of density 
profiles which are built non-parametrically 
from the differentiation of several mass pro- 
files spanning the allowed 68% confidence re- 
gion. In building these density profiles we im- 
pose the physical condition of a negative log- 
arithmic derivative of the mass density at all 
radii. 

In order to determine the 3-dimensional 
number density profile, v (and its confidence 
intervals), up r = 2r u , we use MAPELN. In 
this case we do not exclude late-type spiral 
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Fig. 8. — Mass profiles determined by the 
caustic method. Top panel: the mass profile 
determined using the function (solid line) 
within its 68% c.l. (dotted lines). The dashed 
line shows the mass profile determined using 
Tq = 0.5. Thin solid lines truncated at r = r u 
indicate the 68% c.l.'s on the mass profile de- 
termined by the Jeans analysis (see § 3.1). 
Bottom panel: the combined mass profile de- 
termined via the Jeans and the caustic meth- 
ods. The two profiles at r < r u are obtained 
by integration of the best-fit £- and (3- density 
models (solid and dashed lines, respectively). 
The two profiles at r > r u are obtained via the 
caustic method, by imposing continuity with 
the £ mass profile at r = r u , assuming either 
a variable or a constant (dash-dotted and 
long-dashed lines, respectively). Dotted lines 
indicate the 68% c.l.'s on the mass profile de- 
termined via the Jeans analysis at r < r u , and 
via the caustic analysis at r > r u (continuity 
is imposed at r = r u ). 



members, since we are interested in the to- 
tal number of galaxies in the cluster, irrespec- 
tive of their kinematics. Since the 2dFGRS 
is complete in apparent magnitude, and our 
43 clusters span a significant redshift range 
(0.02-0.13), we need to impose an absolute 
magnitude limit to ensure the same sampling 
of the luminosity functions of all our 43 clus- 
ters. We therefore select the 399 cluster mem- 
bers brighter than the k-corrected (Madg- 
wick et al. 2002) absolute magnitude = 
— 19.3 + 5 log h, which is the magnitude of the 
faintest member of the most distant among 
our 43 clusters, APM 294 (see § 2). This mag- 
nitude is sufficiently bright as to avoid pos- 
sible problems related to the increase in the 
number of dwarf galaxies with clustercentric 
distance (see, e.g. Beijersbergen et al. 2002; 
Durret et al. 2002; Lobo et al. 1997). 

Since the cluster membership assignment 
is problematic beyond r u , we correct the 
MAPELN density profile at larger distances, 
by making use of the optimal caustic. Specif- 
ically, beyond r u we multiply the MAPELN 
density profile by the fraction of cluster mem- 
bers which are contained within the optimal 
caustic in each given radial bin. This fraction 
is 0.8 on average, i.e. ~ 20% of our (pre- 
sumed) cluster members beyond r u could be 
interlopers. 

The mass-to-number density ratio 
function of the clustercentric radial distance is 
shown in Figure 9 (upper panel) . The vertical 
scale in this plot is in normalized units (see 
§ 3.1). The irregular behavior at radii larger 
than r u is due to the use of a non-parametric 
mass profile estimator (the caustic method). 
The bump at r ~ r u is related to the fact that 
we have joined the mass profile determined 
from the Jeans analysis and the mass profile 
determined from the caustic method, with- 
out imposing a continuous derivative. Apart 
from obvious irregularities, one can notice an 
almost monotonous decrease from r ~ 0.3 r u 
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0.0 0.5 1.0 1.5 2.0 

r/r u 

Fig. 9. — Mass-to-number density profiles for 
the ensemble cluster (solid lines) and best- 
fit constants (dashed lines). Top panel: the 
number density profile is determined for the 
sample of all cluster members brighter than 
= —19.3 + 51og/i. Bottom panel: the 
number density profile is determined for the 
sample of early-type cluster members brighter 
than Mbj = —19.3 + 51og/i. The error bars 
account for the uncertainties in both the mass 
and the number density profiles. The vertical 
scale is in normalized units (see § 3.1). 



to the outermost radius. We test the null hy- 
pothesis of a constant mass-to-number density 
profile by the \ 2 method. The null hypothe- 
sis is fully consistent with the data (3.7% c.l.) 
when only the inner r < r u region is consid- 
ered, but it is rejected (99.9% c.l.) over the 
full radial range (0 < r < 2r u ). Had we not 
corrected for the presence of doubtful cluster 
members at radii larger than r u , the decrease 
would have been (marginally) even stronger. 

In Figure 9 (lower panel) we also plot the 
mass-to-number density ratio for the sam- 
ple of 299 early-type, i.e. rj < — 1 (Madg- 
wick et al. 2002), cluster members brighter 
than Mb j = —19.3 + 51og/i. The mass- 
to-number density profile is somewhat flat- 
ter, as expected, since the early-type clus- 
ter galaxies have the steepest number den- 
sity profile among cluster galaxy populations 
(e.g. Dressier 1980; Biviano et al. 2002). In 
this case, the mass-to-number density profile 
is consistent with being constant, not only 
within the central r < r u region (0.2% c.l.), 
but also up to 2r u (74.2% c.l.). 

5. DISCUSSION 

In our analysis we have determined the 
mass and mass-to-number density profiles of 
galaxy clusters out to ~ 2r2oo- We find that 
both cuspy profiles and profiles with a core 
are acceptable. Best-fits to the inner ^ V2oo 
region are found for a /3-profile (see eq. 4) with 
[3 = 0.8 and r c ~ 0.06r2oo 5 or a ^-profile (see 
eq. 3) with £ = 1.4 and a ~ 0.33 r2oo- Ac- 
ceptable fits are also found for a King (1962) 
profile with r c ~ 0.14 r2oo and a NFW profile 
with a ~ 0.18r2oo- Density profiles with very 
large core-radii (r c ^ 0.3r2oo) are ruled out 
(90% c.l.). The density profiles that provide 
the best-fits to the data in the inner ^ r2oo 
region, also provide a reasonable description 
of the cluster gravitational potential at larger 
radii, out to ~ 2r2oo- 
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Let us compare our results with previous 
mass-profile determinations. From the analy- 
sis of the CNOC cluster survey data-set Carl- 
berg et al. (1997a) and vdMOO conclude that 
the average mass profile of galaxy clusters is 
well described by the NFW's analytic form, 
with a scale a ~ 0.25 r2oo, although other 
profiles are equally acceptable. Their best- 
fit NFW scale is slightly larger than ours, but 
not significantly so. Note anyway that the 
CNOC clusters are on average more massive 
than the 43 clusters in our sample, so they are 
expected to be less concentrated and have a 
larger NFW scale. 

Geller et al. (1999) find that the NFW pro- 
file provides a good fit to the Coma mass pro- 
file, while a softened isothermal sphere does 
not. Similarly, Rines et al. (2000) are able to 
rule out a singular isothermal sphere for the 
mass profile of the cluster A576, while a NFW 
profile is a good fit to the observed mass pro- 
file. Both the results of Geller et al. (1999) 
and those of Rines et al. (2000) are based on 
the caustic analysis. 

Similar conclusions about the cluster mass 
profile are reached by Markevitch et al. 
(1999) using X-ray data. Both the cluster 
A2199 and the cluster A496 have mass pro- 
files remarkably well approximated by NFW 
models, and deviate from the isothermal mass 
profile (which is too steep). On the other 
hand, Ettori, De Grandi & Molendi (2002) 
find that about half of 22 nearby clusters of 
galaxies have an X-ray determined mass den- 
sity profile which is better fitted by a King 
(1962) 's rather than a NFW model. In their 
X-ray analyses of 12 galaxy clusters, Durret et 
al. (1994) find that the King profile provides a 
good fit to the mass density distribution, but 
with very small core-radii (~ 15 h^ 1 kpc). 

Useful constraints are derived from the 
gravitational lensing analysis. Most au- 
thors (Clowe 2000; Clowe & Schneider 2001; 
Lombardi et al. 2000; Bautz, Arabadjis, k, 



Garmire 2002) conclude in favor of a NFW 
mass profile and against an isothermal sphere 
mass distribution, except perhaps Sheldon et 
al. (2001). 

Our results are in agreement with the 
claims in support of a NFW mass density 
profile and against both a softened 1 and a sin- 
gular isothermal sphere mass profiles, which 
both increase faster with radius than our ob- 
served cluster mass profile. We stress however 
that our data do not rule out the presence of 
a core in galaxy clusters, although this core 
has to be very small, of the order of the size of 
a galaxy. Interestingly, our formal best-fit for 
a NFW profile is obtained for a ~ 0.18r200; 
which is mid-way between the formal best-fit 
obtained for the CNOC clusters by vdMOO 
and the formal best-fit obtained for the poor 
galaxy systems by Mahdavi et al. (1999). 
Since our clusters are of intermediate mass 
between the CNOC clusters and the poor sys- 
tems of galaxies of Mahdavi et al. (1999), a 
scaling of a with the mass of the system is 
suggested, as expected from numerical simu- 
lations (NFW). 

Our most intriguing result is the significant 
decrease of the mass-to-number density ratio 
at radii greater than ~ 0.3r2oo- The mass per 
galaxy at r2oo < r < 2r2oo is only about half 
the mass per galaxy in the cluster central re- 
gions. Although we do not derive the mass-to- 
light profile (which would require an analysis 
of the cluster luminosity function and it is be- 
yond the scope of this paper), this is unlikely 
to be constant. In fact, for the mass-to-light 
profile to be flat, a strong luminosity segrega- 
tion over the entire radial range out to ~ 2r2oo 
is required. Apart from the segregation of the 
very bright galaxies, which only concerns the 
cluster cores (e.g. Biviano et al. 2002), lumi- 
nosity segregation in clusters is also observed 

1 Geller et al. (1999)'s softened isothermal sphere mass 
profile corresponds to our /3-profile with /3 = 2/3. 
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change of the number ratio of dwarf to gi- 
ant galaxies with radius(e.g. Beijersbergen et 
al. 2002; Durret et al. 2002; Lobo et al. 1997). 
However, the luminosity segregation of dwarf 
galaxies is not relevant in our case, since our 
number density profile is determined only on 
galaxies brighter than = —19.3 + 51og/i, 
which is only ~ 0.3 magnitudes fainter than 
M* in the Schechter luminosity function for 
early- type galaxies (Madgwick et al. 2002). 
We do not find evidence of luminosity segrega- 
tion in our sample of members brighter than 
Mbj = —19.3 + 5 log h, as their absolute mag- 
nitudes and clustercentric distances are un- 
correlated (52% c.L, according to a Spearman 
rank correlation test). 

Our mass-to-number density profile implies 
a mass distribution more concentrated than 
the galaxies. Since the hot X-ray emitting 
gas is less concentrated than the galaxies (e.g. 
Eyles et al. 1991; Durret et al. 1994; Cirimele 
et al. 1997), our result also implies that the 
dark matter distribution is more concentrated 
than the baryonic matter distribution. 

Previous determinations of the mass-to- 
light (or mass-to-number) profiles have pro- 
vided evidences both in favor (Carlberg et al. 
1997a; Cirimele et al. 1997; Mahdavi et al. 
1999; van der Marel et al. 2000; Rines et al. 
2001) and against (Koranyi et al. 1998; Rines 
et al. 2000; Carlberg et al. 2001) a constant 
mass-to- light profile. Koranyi et al. (1998) 
find an increasing mass-to-light profile in the 
poor cluster AWM7, from the center to ~ 0.2 
h~ l Mpc. Such an increase is also seen in 
our data (see Figure 9) but it is not signifi- 
cant. Rines et al. (2000) find that the mass-to- 
light profile decreases with clustercentric ra- 
dius, in consistency with our finding. Carl- 
berg et al. (2001) find an increasing mass- 
to-light profile in galaxy groups. Clearly, a 
general consensus on the relative distribution 
of dark and baryonic matter in galaxy sys- 
tems is still lacking. Part of the discrepancy 



among different authors can arise from the 
use of different photometric bands (Rines et 
al. 2001). Cluster members are redder in the 
cluster center, so that blue-band selected clus- 
ter samples (like ours) should be character- 
ized by flatter number density profiles, and 
hence steeper mass-to- number density profiles, 
as compared to samples selected in the near- 
infrared. Our mass-to-number density profile 
indeed becomes flatter, and consistent with 
being constant, when only early-type (and 
hence redder) cluster members are considered 
(see Fig. 9). 

Does our result invalidates the use of sim- 
ple mass estimators based on the light traces 
mass hypothesis? Probably not. In fact, our 
mass-to-number density profile does not sig- 
nificantly deviate from a constant when we 
restrict the analysis to the virialized region 
(r < r2oo), or when only the early- type galax- 
ies are considered. 

6. SUMMARY 

We use the June 2001 public release of the 
2dFGRS to construct a sample of 4947 galax- 
ies with redshifts in the region of 43 clusters, 
reaching out to ~ 2r2oo- 1345 cluster mem- 
bers are selected with the method of Fadda 
et al. (1996). We build an ensemble cluster 
profile from this sample, using r2oo and o p 
of each cluster to normalize the clustercentric 
radii and velocities of all galaxies. 

We determine the mass profile of this en- 
semble cluster by two independent methods, 
namely those of vdMOO and D99. The method 
of vdMOO is based on the solution of the Jeans 
equation, and can only be applied in condi- 
tions of dynamical equilibrium. We therefore 
apply it on the subsample of cluster mem- 
bers within the virial radius, r2oo- We ex- 
clude late-type spirals from this sample since 
these galaxies have been shown in the litera- 
ture to be characterized by anisotropic orbits 
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(e.g. Mohr 1996; Biviano et al. 1997). The 
method of D99 can usefully be applied to the 
whole sample of galaxies in the cluster region, 
since it does not require dynamical equilib- 
rium. We use it to constrain the cluster mass 
profile beyond the virial radius, and out to the 
limit we imposed on our sample, ~ 2r2oo- 

We find very good agreement between the 
mass profiles determined using the two meth- 
ods. The mass profile from the center to 
~ 2r2oo is remarkably well fitted by a £-model 
(eq. 3) with a ~ 0.33 r2oo and £ = 1.4, but 
many other models provide acceptable fits, in- 
cluding NFW models and models with a small 
core (r c ^ 0.2r2oo)- We find an uncertainty on 
the total mass of « ±8% at the virial radius, 
and ±15% at two virial radii. 

We determine the mass-to-number density 
profile of the ensemble cluster. It is fully 
consistent with a constant from the center to 
~ r 200 ; but it significantly decreases at larger 
radii. The decreasing trend is however not 
significant when the number density profile is 
computed on the early-type galaxies only. 

This work is based on the 100k Data Re- 
lease of the Two Degree Field Galaxy Redshift 
Survey, kindly made available to the astro- 
nomical community by the 2dFGRS Team. 
We acknowledge useful discussion with An- 
tonaldo Diaferio, Peter Katgert and Alain 
Mazure. We thank an anonymous referee for 
her/his constructive comments. We thank 
Tim Beers and David Merritt for provid- 
ing us with copies of their FORTRAN codes 
ROSTAT and, respectively, MAPELN. This 
work was partially supported by the Italian 
Space Agency (ASI), and by the Italian Min- 
istry of Education, University, and Research 
(MIUR grant COFIN2001028932 "Clusters 
and groups of galaxies, the interplay of dark 
and baryonic matter"). 
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Table 1 — Continued 
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